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A phase-field model that takes into account the bending energy of fluid vesicles is presented. The 
Canham-Helfrich model is derived in the sharp-interface limit. A dynamic equation for the phase- 
field has been solved numerically to find stationary shapes of vesicles with different topologies and 
the dynamic evolution towards them. The results are in agreement with those found by minimization 
of the Canham-Helfrich free energy. This fact shows that our phase-field model could be applied to 
more complex problems of instabilities. 
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I. INTRODUCTION 

Biological membranes have been of wide interest to 
biologists, chemists and physicists for many years (cfr. 
for a historical review). They are the frontier which 
defines cells and most of their internal organelles, and 
they act as a boundary between internal cellular orga- 
nization and its surrounding medium. Membranes are 
composed of several kinds of lipids (phospholipids, choles- 
terol and glycolipids), which are self-assembled in a fluid 
bilayer, and by membrane proteins which are anchored 
on it JJ. From the molecular point of view, biomem- 
branes are extremely complex. However, there seems to 
be a universal construction principle common to all ac- 
tual membranes, which is the presence of a fluid lipid 
bilayer through which proteins can diffuse. Vesicles are 
model closed membranes consisting of one or several dif- 
ferent kinds of lipids 0, 0| . They have therefore been 
studied to get an idea of the main physical properties of 
actual biomembranes 0, such as red blood cells's jj. 

Recently several experimental results have been re- 
ported on dynamic instabilities in membranes, such as 
pearling (I^ ^ , budding and tubulation ^ . In these 
experiments shape instabilities are induced by the inser- 
tion of a certain concentration (locally or globally) of an 
amphiphilic polymer (which mimicks the proteins within 
the biomembranes) in the outer leaflet of the bilayer [lOl| . 

Dynamic models of vesicles in shear flow ^J, ll^ 
IT^ . domain growth in multicomponent membranes and 
bud formation 

El III III 111 , and pearhng [3 

have 

been studied theoretically. The effect of proteins included 
in lipidic membranes is also the subject of theoretical 
studies Umilll. 

Phase-field models (or diffuse interface models) can 
be thought as mathematical tools to study complex 
interfacial problems, such as free boundary problems 
p^. Phase- field models are mesoscopic models of the 
Ginzburg-Landau type, which disregard microscopic de- 
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tails. Such models have been widely used before in dif- 
ferent interfacial problems such as solidification and the 
Saffman- Taylor problem and roughening ||2J]. Most 
of these phase-field models describe the effect of surface 
tension, but do not deal with bending energies. 

Biben and Misbah studied the tumbling transition of 
a vesicle under shear flow ^3 using a phase-field-like 
model. Their model couples the bending energy with 
the velocity field as a force in the two dimensional hy- 
drodynamic equation (in the Stokes limit). Du et al. 
used a phase-field model |2l |2^ based on the free ener- 
gies for surfactant films found in the works of Gompper 
and Schick . They applied their model to min- 

imize the free energy to find stationary shapes, using a 
standard global Lagrange formalism to deal with both 
global area and volume constraints. The minimization 
was performed by a gradient fiow equation, through a 
non-conserved purely relaxational dynamics. 

In this paper we have derived a phase-field model for 
the bending energy of fluid vesicles. The membrane 
is considered as a mathematical interface between two 
phases, the inner fluid and the outer fluid. In this kind of 
models there is no need to track the interface during the 
dynamic evolution, which is one of the main problems in 
membrane dynamics Our equations are continuous 
in the whole domain, and the interface is located by the 
level-set of the phase-field, i.e. the region of rapid varia- 
tion of the phase-field. The free energy functional asso- 
ciated with this model reduces to the Canham-Helfrich 
bending energy of the lipid bilayer |2l in the so- 
called sharp interface limit, when interface width goes to 
zero. In addition, phase-field models are dynamic mod- 
els, so we are capable with our model to study dynamic 
properties of vesicles, such as relaxation towards station- 
ary shapes. The fact that we find the correct stationary 
shapes shows that our free energy functional deals cor- 
rectly with bending energies. 

This paper is organized as follows: In Sec. ^] we pro- 
pose a phase-field model to take into account the bending 
energy of fluid vesicles, and we prove that it is equiva- 
lent to the Canham-Helfrich model. An effective free 
energy functional is written down in the presence of geo- 
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metric constraints, and the dynamic equation is worked 
out. Details of the numerical integration are presented 
in Sec. IIIII The stationary shapes for spherical and non- 
spherical topologies are presented in Sec. IIVI together 
with some features of dynamic evolution. A shape di- 
agram is presented in Sec. for vesicles with spherical 
topology. Discussion of the possible applications and ex- 
tensions of this model, and main conclusions are found 
at the end of the paper in Sec. IVII 

II. MODEL 

A. Canham— Helfrich free energy 

A mathematical model describing the geometric prop- 
erties of lipidic membranes was introduced by Can- 
ham '2^ and Helfrich 30] to describe interfaces governed 
by bending energy. The membrane is considered to be 
a two-dimensional surface embedded in Euclidean three- 
dimensional space. Such a surface has two curvature in- 
variants: the mean and the Gaussian curvatures. The so- 
called Canham-Helfrich local bending energy is thus an 
expansion in terms of these invariants. For membranes 
with a possible assymmetry between the two layers, the 
spontaneous curvature model Q reads 

/C-H, sc = I (2i? - Co)' + KG^, (1) 

where k and are two elastic constants: the bending 
rigidity, and the Gaussian bending rigidity, respectively; 
and Co is the spontaneous curvature. Due to the Gauss- 
Bonnet theorem, the Gaussian curvature term (the last 
term in Eq. QJ) integrated over a closed surface is a topo- 
logical invariant. Since we are not concerned with study- 
ing topological changes here, this term will be a constant 
factor in the total free energy, so it does not need to be 
considered. Therefore the bending energy reduces (after 
integration over the whole vesicle membrane surface) to 

i^C-H, sc = |^(2H-co)'ds, (2) 

where T is the membrane surface. The special case where 
Co = in Eq. (jSJ is the minimal model 

Fc-H^^J^i2Hfds, (3) 

where no asymmetry between the two layers of the lipidic 
membrane is considered. 

B. Phase— field implementation 

1. Minimal Model 

In order to obtain our phase-field model, we will start 
with an ansatz for the free energy functional, which will 



be shown to be equivalent, under certain limits, to the 
Canham-Helfrich Hamiltonian Eq. (|3J), with no sponta- 
neous curvature. This free energy functional is 

m = / da;^ (4) 

Jn 

where the function is the so-called phase-field, and 

$[0(a;)] = -(/) + 0^-e2v20(a;), (5) 

where e is a small parameter related to the interface 
width. Note that the free energy density functional in 
Eq. is nothing else other than the square of the 

chemical potential (the functional derivative of the free 
energy) associated with the Cahn-Hilliard problem (3ll |. 

The minimum of the free energy Eq. Q), with no 
constraints, is obtained by setting Eq. (|SJ) equal to 
0. In one dimension, this leads to the tanh-like solu- 
tion 4){x) — tanh (y^-^j^^ J given the boundary conditions 

(/)(±cxd) = ±1. The boundary conditions in three dimen- 
sions are that the phase-field at infinity is = — 1, which 
is the value for the stable phase of the outside bulk. 

Phase fields are regular functions, so they can be writ- 
ten in terms of the signed distance to the interface, the 
set (j) = Q, in units of e, 

cj>{x) = f{D{x)), (6) 

where D{x) = d{x)/e, and d{x) is the signed distance of 
the point x to the surface F. We should plug this into 
Eq. However, before doing this, we will explicitly 

work out the Laplacian term 

V^* = v(vf(D{x))) = V (f'(D(x))^^^ 

We now consider that our surface is regular enough to 
let the distance be just the position vector (plus some 
constant vector, which will be irrelevant when differenti- 
ated). Thus, 

Vd{x) = h, (8) 

is a unit vector normal to the interface. If we insert 
Eq. ijSJl into Eq. Q, we obtain 

= ^.f"{D{x)) + -J'{D{x))V^d{x). (9) 

Thus, taking the result from Eq. we see that Eq. @ 
reduces to 

Hf)^ I [{!" ^[f~l)l)+ef{D[x))V^d{x)]'dx. 

(10) 

This is the free energy functional that we would like to 
minimize dynamically under certain constraints. 
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In the sharp-interface hmit (e — > 0) and, in order to 
minimize the free energy functional, Eq. I|1U|) . we obtain 

r{Dix)) {f{Dix)) 1) f{D{x)) ^ 0. (11) 



Therefore we obtain the interfacial solution 

'd{xy 



f{D{x)) = tanh 



(12) 



where we can see that e\/2 is the width of the interface. 
Thus, Eq. (Uni reduces to 



Hf)= f \eriDix))V'd{x)\^dx. 
Jn 



From Eq. H12|l . we can compute 
/' iD{x)) = i=sech2 



d{x) 

77f 



(13) 



(14) 



Using the fact that 
lim 



sech 



d{x) 
eV2 



S{d{x)), (15) 



we can write the integral over the volume Q in Eq. (|13|l . 
as an integral over the surface F 



^(rf) = ^ ^3 f \V^dix)\^ds. 

3 Jr 

Using the fact that (cfr. e.g. |3^. l3^) 



H = -V'd{x), 



(16) 



(17) 



we can see that the free energy Q is, up to a known pro- 
portionality constant, the same as the Canham-Helfrich 
minimal model Eq. We may thus employ the phase- 
field dependent functional to implement the bending en- 
ergy in a vesicle membrane. Since the proportionality 
constant is positive, it just redefines the energy scale. It 
is thus exactly equivalent to minimizing our free energy 
functional instead of the curvature dependent one. 



2. Spontaneous Curvature 

It is possible to extend the results of the last section in 
order to include a spontaneous curvature, as in Eq. 
The new ansatz for the free energy functional will be the 
same as in Eq. Q), but with an additional term in the 
free energy density 



- / 

In 



dx, 



where 



ix)]^^^{x)]+eCoix)il-^^), 



(18) 



(19) 



where Co{x) is, as we show below, related to the spon- 
taneous curvature which, in principle, can be position- 
dependent (or even 0-dependent). The proof that Eg. 1181 
is equivalent to the SC model Eq. ^ is analogous to that 
shown in Sec. Ill Bl Therefore, Eq. (|10l) will read now as 

•^.c(/) = / [{r-{f-i)f) + e{rv'dix) 

Jn 

-C^{l~f))]'dx. (20) 
The order e term can be rewritten as 
/' V^d{x) - Co(l - f) = (V^d{x) - V2C^ /', (21) 



where we used the fact that \ — p = V2/', which comes 
directly from Eq. H12|l . We finally reach 



ds. 



(22) 



which, identifying ^/2Co = cq, is equivalent (again up to 
a constant) to the spontaneous curvature model Eq. 1^)1. 



C. Geometrical constraints 

Vesicle shapes are subject to certain geometrical con- 
straints. Their enclosed volume and their surface area 
should remain constant. We will need to implement these 
constraints in our phase-field model. 



1. Local surface area 

Lipidic membranes are in a fluid phase at physiolog- 
ical temperatures. The solubility of membrane lipids is 
extremely low, which implies no relevant exchange of ma- 
terial between the membrane and the surrounding media. 
In addition, the membrane can be considered to be lo- 
cally incompressible. These two facts provide us with a 
constraint for the local area of the vesicle, which remains 
fixed. 

We have implemented this constraint in our phase- field 
model via a Lagrange multiplier function coupled with 
the surface area in the free energy functional. Thus, we 
define an effective free energy functional 



^cff[0] = ^[0] + / (j{x)a[(t>\dx, 
Jn 



(23) 



where T[4>\ is given by Eq. Q, cr is a Lagrange multiplier 
(interpreted as a surface tension), and a[(/)] is the local 
surface area functional, 



a{x) — 6 {d{x)) 



(24) 



which we rewrite in terms of the parameter e using the 
representation of the delta, Eq. (|15|l . 



2V2 



(25) 
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where we have used the fact that e'V(j){x) = f'{D{x)) h, 
and, from Eq. that {f'{D{x)))^ - S{d{x)). Using 

Eq. (|24|l . i.e. Eq. 1)25(1 in the sharp-interface hmit, the 
last expression integrated over the whole domain Q is 
equivalent to the surface area of the vesicle, 



where we have defined a as 



a \(h] dx 



ds. 



(26) 



2. Enclosed volume 

The other geometrical property of the vesicle that re- 
quires a constraint is its enclosed volume. Biological 
membranes are permeable to water, but not to, e.; 
large ions (on the time scales we are interested in) 
This means that any transfer of water through the mem- 
brane would create an osmotic pressure which cannot be 
counterbalanced by the relatively much weaker bending 
energy, 30]. Therefore, the concentration of osmotically 
active molecules fixes the inner volume of the vesicle. The 
usual way to implement this condition in the free energy 
is, as done before with the surface area constraint, to in- 
troduce a Lagrange multiplier coupled with the volume 
term, which ensures its conservation. 

Thus, dynamics are introduced such that the vesicle 
volume is conserved through time. Model-B-like dynam- 
ics is therefore used 



dt \ 6(j} 



(27) 



This dynamic equation ensures that (l)(x)dx is con- 
stant in time. This integral is equal to the difference of 
the inner and outer volumes (<f> takes its stable values -1-1 
and —1 inside and outside the vesicle, respectively), as 
e 0. As the sum of the inner and outer volumes is the 
volume of the space £7 (which is constant), then we can 
write the inner volume as 



v{n) 



{x)dx 



(28) 



where it can be seen that model-B dynamics ensures con- 
servation of the inner volume through the dynamic evo- 
lution. 



D. Dynamic equation 

The dynamic relaxation towards free energy minima is 
achieved in our model by conserved relaxation dynamics, 
Eq. (|27|l . We could have introduced some homogeneous 
mobility, but this would have just changed the time scale. 
Performing the functional derivative in Eq. H27f) gives the 
dynamic equation 

^ = 2V^[{^ct>''-l)^^-^V^m+^a{x)V^ 



+ e2Va(a;)- V0}, 



a{x) = ^cr(a?). 



(30) 



Using this kind of dynamics, local conservation of the 
inner volume of the vesicle is achieved in a natural way, 
unlike Ref. '26"! which uses a purely relaxational dynamics 
with no direct conservation of the inner volume of the 
vesicle. 



III. NUMERICAL INTEGRATION 

The partial differential equation (|29|l is highly non- 
linear (notice, e.g., the coupling between the field 0^ and 
the functional $[</)]). We thus need to use a numerical 
method to solve this equation in order to study dynamic 
relaxation towards stationary vesicle shapes. 

The discretization algorithms used are second-order fi- 
nite differences for the spatial dependence, and an Euler 
scheme for the time dependence '34]. 

Our effective free energy functional H23|) explicitly con- 
tains a Lagrange multiplier. Therefore, we need to know 
the time evolution of the Lagrange multiplier. To do this, 
we have used a first order Lagrangian method to study 
how the multiplier evolves to its stationary value [ssj ] 



{x) = a'^ix) + a (aict^^ix)] ~ ao{x)) , (31) 



(29) 



where a > is the stepsize, k is the discretized time, and 
ao{x) is the fixed local surface area. Since we are not 
interested in the actual dynamics of the multiplier, our 
choice is justified because it does not change the dynam- 
ics of the phase-field, but it just keeps the surface area 
of the vesicle constant during the time evolution without 
altering the dynamics. 

We have performed simulations on lattices of differ- 
ent sizes and equivalent shapes and evolutions were ob- 
tained. As the discretization of the differential equation 
(|29|l was done with a consistent finite difference method 
(standard second-order finite differences), and the time 
step was chosen following the Courant-Friedrichs-Lewy 
stability criterion. At < \k\ Ax, where k is some con- 
stant, we can thus assume that the algorithms used are 
convergent [s^. In addition, during the time evolution, 
we checked the value of the free energy evolution in time 
to see how it relaxes to a stationary value in a monotoni- 
cally decreasing way. The values of the inner volume and 
the surface area were also computed during the evolution 
and it can be seen that the volume remains constant (up 
to the numerical precision) during all the process, and 
similarly with the surface area (the value of the Lagrange 
multiplier converges rapidly to the stationary solution). 



IV. RESULTS 

In our model there seem to be several free parameters 
(e, ao{x), Vo). However, e is a small parameter (the model 
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is shown to be robust under variations of this parameter), 
which will be set, in what follows, to be equal to the mesh 
size. 

In addition, the term proportional to ^a{x) in the dy- 
namic equation (the last term in Eq. H29|l) is shown to 
be small, and the Lagrange multiplier, a, can be consid- 
ered homogeneous. To show this, we numerically com- 
pute these variations for an initial profile of a{x), and 
see how they rapidly relaxe to a nearly constant func- 
tion (i.e the time scale of the relaxation of the effective 
surface tension is smaller than the time scale related to 
the shape change). Moreover, a{x) appears as an effec- 
tive surface tension which prevents the surface area from 
changing. Anyway, its value in membranes is very small 
compared with other energy scales in the system (e.g. 
bending rigidity) |3(i| . Therefore, its variations are also 
small. 

Finally, scale invariance causes that the ratio between 
the constrained total volume and the total surface area 
is the only relevant parameter in the model (for a fixed 
topology). Thus, we define a dimensionless volume v as 
the ratio between the actual volume and the volume of a 
sphere with the same area. 



V 



where 



1/2 



(32) 



(33) 



We will thus look for stationary shapes for different fixed 
topologies: spherical (Euler characteristic equal to 2) and 
non-spherical (e.g. genus-1 toroidal topologies with Eu- 
ler characteristic equal to 0). We will focus mainly on 
spherical topology, in order to discuss the model and the 
results obtained. 



A. Spherical topology 

For this topology, the three qualitatively different sta- 
tionary shapes for the minimal model are found, in agree- 
ment with [33 • These are, in order of decreasing reduced 
volume, the prolate and oblate ellipsoids and the stoma- 
tocyte (see Fig. Pl. 

Shape evolution is done starting from an arbitrary ini- 
tial shape. Since the dynamic equation has no external 
noise (just numerical noise), we start from different initial 
shapes corresponding to each value of the reduced vol- 
ume. The initial condition for the phase-field is a sharp 
distribution of (/> = -I- 1 and (j) = —1. There is thus a tran- 
sient period in the first few steps of the evolution, where 
the diffuse interface is created and a tanh-like profile is 
obtained, which remains during the subsequent evolu- 
tion. This transient is needed to calculate the surface 
area using Eq. (|25|l . since a gradient in is required. 




(a) V = 0.69 





(b) V = 0.60 



(c) V = 0.43 



FIG. 1: Stationary shapes for the minimal model, as station- 
ary states of the dynamic evolution of certain initial condi- 
tions under Eq. 12911 . (a) Prolate, (b) oblate, and (c) stoma- 
tocyte are shown. 



We have solved Eq. (|29|l numerically on a three- 
dimensional lattice. The possibility of finding non- 
axisymmetric shapes then arises. In Fig. |21 four snap- 
shots of the shape evolution towards a prolate ellipsoid 
with V = 0.69 are shown. This is a stable shape, since the 
actual transition between prelates and oblates happens at 
a value of the reduced volume that is lower than this value 
{vd — 0.65) [3^ . The dynamic evolution towards this ax- 
isymmetric prolate is done with a non-axisymmetric dy- 
namics, and no axis of symmetry has been supposed. The 
initial shape (Fig.|2Ia)) is a non-axisymmetric box, which 
dynamically evolves towards an axisymmetric shape. 

We have also studied the behavior of the dynamic equa- 
tion in the axisymmetric case, where we can discretize 
Eq. H29|l on a two-dimensional lattice. Fig. O shows 
the time evolution to eventually reach a stomatocyte-like 
shape. Plots show the value of the phase-field on a grey- 
level scale, where black represents 4> = +I7 or the inner 
volume of the vesicle; and white corresponds to (j) = ~lj 
or the outer volume of the vesicle. 

We have also found the level-set of these grey-level 
scale plots, to track the position of the vesicle membrane. 
In Fig. 2] a continuous fit of this contour is plotted and 
the evolution towards a discocytc-like shape is shown. 
We can see that the reduced volume corresponding to 
this figure is v ~ 0.5, which is conserved during evolu- 
tion (so they are the surface area and the volume sepa- 
rately). This shape is known to be metastable. However, 
this shape is obtained from the initial ellipsoid because 
the actual stable shape, corresponding to that value of 
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(a) 10"' A< (b) 0.5 X lO^^ At 




(c) 2 X 10^ At (d) 8 X IQS At 




(c) 3.3 X 10** At; V ~ 0.492 (d) 7.8 X 10** At; v ~ 0.496 



FIG. 4; Snapshots of the time evolution towards a discocyte- 
like shape. Figures shown here are for the axisymmetric case. 
The dot-dashed hue indicates the axis of symmetry. Curves 
are continuous fits of the level-set of the phase-field on a 80 x 
40 lattice. The time step is set to At = 10~*. 



FIG. 2: Shape evolution for a vesicle with v = 0.69, which 
eventually reaches a prolate shape, the integration was per- 
formed on a 3-dimensional 50 x 50 x 50 lattice with time step 
At = 10-*. No axis of symmetry is supposed here, and the 
initial shape is a non- axisymmetric 40 x 10 x 10 box. 



states where the value of the energy remains constant. 



B. Non-spherical topologies 



1 




> 


(a) 10^ At 


(b) 0.5 x lO'^ At 









(c) lO'^ At (d) 8 x lO'^ At 



FIG. 3: Shape evolution for a vesicle with v — 0.43, which 
eventually reaches a stomatocyte shape. A two-dimensional 
section is shown, where an axis of symmetry exists, located on 
the lower side of each snapshot. Integration was performed 
on an axisymmetric 60 x 30 lattice and the time step was 
At = lO"''. 



the reduced volume, is far away in the shape landscape 
from our initial choice. Note also that the dynamic equa- 
tion H29|) is such as the free energy is a monotonically 
decreasing function, which reaches metastable or stable 



We have also studied non-spherical topologies, such as 
the genus-1 toroidal topology [3^ (see Fig. [Sfa) and 
(b)). Circular tori are found for large values of the re- 
duced volume V, and sickle-shaped tori for small values 
of V. Discoids can be found for intermediate values of w, 
although they are not stable shapes, and they will even- 
tually fall to stable ones. In addition, spherical shells are 
also found (see Fig. Efc)). They have a different topol- 
ogy with an Euler characteristic % = 4. These shapes 
can be thought of as a limit case of a sickle-shaped torus, 
when the outer radius vanishes thus changing the global 
topology of the shape (Fig. |SJ). 

The shapes shown in Fig.|31may be found in different 
ways. First, it is possible to take an initial shape of a 
given topology, and let it relax to a stationary shape of 
that same topology. Within this relaxation there is no 
topological change, and a non-spherical stationary shape 
may be found. Second, in some cases, dynamic evolution 
within a given topology leads to shapes close to a topo- 
logical transition. Then, because of the natural way of 
dealing with topological transitions of phase-field mod- 
els, topology may change. We show an example of how 
we get the shape in Fig. Hfc) in this way in Fig. |H1 In 
this Figure there is a topological transition between steps 
Efd) andinie). The poles of the torus get closer to the axis 
of symmetry and, eventually, fuse to get spherical-shell- 
like topology. The actual dynamics of this transition is 
not explained by our model, neither by energetics of the 
Canham-Helfrich model, since the process of membrane 
fusion is far from being well-understood. When mem- 
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branes are fused or broken, the Helfrich approach is not 
vahd, since the microscopic details of the bilayer become 
then relevant (namely, the radii of curvatures involved are 
of the order of the membrane width). In addition, when 
a topological change occurs, there is an energy change 
due to the Gaussian curvature term in Eq. (jT]). 

In any case, once the shape has changed its topology 
the evolution leads the shape to a stationary state given 
that new topology (we can think of the initial shape being 
the shape after the topological change). 




(c) V = 0.27 



FIG. 5: Stationary shapes for non-spherical topologies, as 
stationary states of the dynamic evolution of certain initial 
conditions under Eq. 1291 (a) Sickle-shaped torus of genus- 
1 toroidal topology, (b) circular toroid (Clifford torus js^) 
and genus- 1 toroidal topology, and (c) spherical shell with an 
Euler characteristic x = 4 are shown. 



V. SHAPE DIAGRAM 

The results shown in the last section are collected in 
a shape diagram where the bending energy is plotted 
against the reduced volume (see Fig. d) . The curvature 
energy of the shapes obtained as stationary states of the 
dynamic evolution under Eq. H29|l . is calculated in the 
following way. Firstly, the interface is located as the 
level-set, (/> = 0, of the phase-field. An interpolation is 
then performed over the discrete data in order to obtain 
a continuous function descri bing the membrane. Using 
surface differential geometry |32L 13^. the curvature ten- 
sor on the surface defined by the interpolating function 
rotated about the axis of symmetry is worked out. The 















(a) 10^ At 


(b) 10'^ At 








w ^ 


(c) 5 X 10'' At 


(d) 6.8 X lO'^ At 


n 




(e) 7 X 10'' At 


(f) 8 X 10^ At 



FIG. 6: Shape evolution for a vesicle with v = 0.27, which 
eventually reaches a spherical shell shape. A two-dimensional 
section is shown, where an axis of symmetry exists, located on 
the lower side of each snapshot. Integration was performed on 
an axisymmetric 50 x 25 lattice and the time step was At = 
10~^. The initial shape (a) has a genus- 1 toroidal topology. 
It evolves dynamically untill it changes its topology towards 
an spherical shell. 



trace of this tensor, related with the mean curvature, is 
then calculated. Integrating over the surface eventually 
gives the bending energy. 

Once these pairs {v, F} are obtained, a comparison 
with known results is presented. Thus, in Fig. [7| we 
plot the lines corresponding to the minimization of the 
Canham- Helfrich free energy, Eq. with fixed vol- 
ume and surface area js^l- There are three branches 
of different shapes (stomatocytes, oblates and prolates), 
which intersect at certain values of the reduced volume, 
V, where a change in the stability of the shape occurs. 
We see that the results obtained for stationary shapes 
with the phase-field Eq. |0J arc in good agreement with 
this. 



VI. DISCUSSION AND CONCLUSIONS 

The model could be extended to include hydrody- 
namic effects. Since in this paper we give a derivation 
of the Canham-Helfrich free energy within a phase-field 
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FIG. 7: Shape diagram for the minimal model with spher- 
ical topology. Lines correspond to the minimization of the 
Canham-Helfrich energy (C-H) (37| . and symbols to the re- 
sults of the phase-field model (P-F). 



approach, a hydrodynamic equation for the velocity field 
(e.g. the Navier-Stokes equation) could be introduced, 
with a force acting on the membrane due to the bending 
elasticity. In addition, an advective term should be in- 
cluded in the dynamic equation for the phase- field 'l^ . 
However, our future aim is to study the shape instabilities 
seen in There, the relevant characteristic time scale 

that needs to be studied is that associated with the relax- 
ation of the curvature, which turns out to be related to 
the diffusion coefficient of the polymer in the membrane, 
and not directly to the membrane viscosity. Therefore, a 
dynamic model which couples the polymer concentration 
in the membrane with the local spontaneous curvature 
would be needed. 

The existence of a differential flow between monolay- 
ers (i.e. sliding between both monolayers) due to the 
viscous drag between the non-polar tails of the hydro- 
carbon chains of the lipids is a microscopic effect which 
may induce a velocity discontinuity at the membrane. 
Thus, it would be an important factor to keep in mind 
when introducing the hydrodynamic field [i^, |41| . When 
dealing with shape instabilities due to the anchorage of 



amphiphilic molecules 0, 0, curvature can be induced 
by two mechanisms, both of which are theoretically ex- 
plained by a model: the area difference elasticity (ADE) 
model, and the spontaneous curvature (SC) model jjj. 
ADE relaxes via the sliding between monolayers, and SC 
via diffusion of the polymer in the membrane. From ex- 
periments, it is possible to distinguish between the ADE 
and the SC mechanisms by measuring the relaxation 
times (or the diffusion constants). In the experiments 
done by Tsafrir et al. 0, evidence that both ADE and 
SC mechanisms influence pearling instability was pro- 
vided. However, it was assumed there that pearling was 
a result of the SC mechanism, i.e. that the dynamic 
instability was not affected by the sliding between the 
monolayers of the membrane. 

As a conclusion, a phase-field model for bending energy 
of fluid vesicles, and not surface tension as in usual phase- 
field models, has been derived. It has been shown to be 
equivalent in the sharp-interface limit to the Canham- 
Helfrich free energy for closed membranes. Various geo- 
metric constraints are also implemented in the phase-field 
approach. Dynamic equations for the phase-field and 
for the Lagrange multiplier related to the surface area 
conservation have been worked out explicitly. These are 
highly non-linear partial differential equations and thus 
numerical treatment was found to be necessary. Station- 
ary shapes of vesicles with different topologies are found, 
in agreement with those obtained by minimization of the 
Canham-Helfrich free energy In addition, they are 
found dynamically, from arbitrary initial shapes. More- 
over, a shape diagram for spherical topology vesicles is 
presented. These facts show that our phase- field model 
is a good description of the bending elasticity of mem- 
branes, and that it could be used within a generalized dy- 
namic framework where hydrodynamics or polymer dif- 
fusion on the membrane could be introduced. 
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